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We study decoherence of a quantum dot charge qubit due to couphng to piezoelectric acoustic phonons in the 
Bom-Markov approximation. After including appropriate form factors, we find that phonon decoherence rates 
are one to two orders of magnitude weaker than was previously predicted. We calculate the dependence of the 
Q-factor on lattice temperature, quantum dot size, and interdot coupling. Our results suggest that mechanisms 
other than phonon decoherence play a more significant role in current experimental setups. 
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I. INTRODUCTION 

Since the discovery that quantum algorithms can solve cer- 
tain computational problems much more efficiently than clas- 
sical ones,' attention has been devoted to the physical im- 
plementation of quantum computation. Among the many 
proposals, there are those based on the electron spin^-^^ or 
chargei^i^iiS in laterally confined quantum dots, which may 
have great potential for scalability and integration within cur- 
rent technologies. 

Single qubit operations involving the spin of an electron 
in a quantum dot will likely require precise engineering of 
the underlying material or control over local magnetic fields 
both have yet to be achieved in practice. In contrast, single 
qubit operations involving charge in a double quantum dot 
(DQD)'" are already within experimental reachiiiii^ They can 
be performed either by sending electrical pulses to modulate 
the potential barrier between the dots (tunnel pulsing)* ** or 
by changing the relative position of the energy levels (bias 
pulsing)." In both cases one acts on the overlap between the 
electronic wave functions of the dots. This permits direct con- 
trol over the two low-energy charge states of the system - the 
basis states |1) and |2) of a qubit: Calling A^i (iV2) the num- 
ber of excess electrons in the left (right) dot, we have that 
|1) = (1,0) and|2) = (0,1). 

The proposed DQD charge qubit relies on having two lat- 
eral quantum dots tuned to the (1,0)^ (0, 1) transition line of 
the Coulomb blockade stability diagram (see Fig.Q. Along 
this line, an electron can move between the dots with no charg- 
ing energy cost. An advantage of this system is that the Hilbert 
space is two-dimensional, even at moderate temperatures, 
since single-particle excitations do not alter the charge con- 
figuration. Leakage from the computational space involves 
energies of order the charging energy which is quite large in 
practice 1 meV ^ 10 K). In the case of tunnel pulsing, 
working adiabatically - such that the inverse of the switching 
time is much less than the charging energy - assures mini- 
mal leakage. The large charging energy implies that pulses as 
short as tens to hundreds of picoseconds would be well within 
the adiabatic regime. However, the drawback of using charge 
to build qubits is the high decoherence rates when compared 
to spin. Since for any successful qubit one must be able to 



perform single- and double-qubit operations much faster than 
the decoherence time, a quantitative understanding of deco- 
herence mechanisms in a DQD is essential. 

In this work, we carry out an analysis of phonon decoher- 
ence in a DQD charge qubit. During qubit operations, the 
electron charge movement induces phonon creation and anni- 
hilation, thus leading to energy relaxation and decoherence. 
In order to quantify these effects, we follow the time depen- 
dence of the system's reduced density matrix, after tracing out 
the phonon bath, using the Redfield formalism in the Born and 
Markov approximations. '■'•ii 

Our results show that decoherence rates for this situation 
are one to two orders of magnitude weaker than previously 
estimated. The discrepancy arises mainly due to the use of 
different spectral functions. Our model incorporates realistic 
geometric features which were lacking in previous calcula- 
tions. When compared to recent experimental results, our cal- 
culations indicate that phonons are likely not the main source 
of decoherence in current DQD setups. 

The paper is organized as follows. In Sec.|llJ we introduce 
the model used to describe the DQD, discuss the coupling to 
phonons, and establish the Markov formulation used to solve 
for the reduced density matrix. In Sec. [Ill] we study decoher- 
ence in a single-qubit operation, while in Sec. lIVI we simulate 



{0,1} \ 


(1,1) \^ ^ 


A . 

9 






(1,0) ^ 



V, 



FIG. 1 : Schematic Coulomb blockade stability diagram for a double 
quantum dot system at zero bias4- {Ni,N2) denotes the number 
of excess electrons in the dots for given values of the gate voltages 
Vi and V2. The solid lines indicate transitions in the total charge, 
while the dotted lines indicate transitions where charge only moves 
between dots. The point A marks the qubit working point. 
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the bias pulsing experiment of Ref.Tl', Finally, in SecMwe 
present our conclusions. 



II. MODEL SYSTEM 

We begin by assuming that the DQD is isolated from the 
leads. The DQD and the phonon bath combined can then be 
described by the total Hamiltoniani 



H = Hs + Hb+Hsb, 



(1) 



where Hs and Hb are individual DQD and phonon Hamil- 
tonians, respectively, and Hsb is the electron-phonon inter- 
action. We assume that gate voltages are tuned to bring the 
system near the degeneracy point A (Fig. where a single 
electron may move between the two dots with little charging 
energy cost. To simplify the presentation, only one quantum 
level on each dot is included; En2) denotes the energy of 
an excess electron on the left (right) QD (possibly including 
some charging energy). Likewise, spin effects are neglected. 
Thus, in the basis {|1), |2)}, the DQD Hamiltonian reads 



s{t) 



(2) 



where cr^.j^ are Pauli matrices, e{t) — Ei — E2 is the energy 
level difference, and v{t) is the tunneling amplitude connect- 
ing the dots. Notice that both e and v may be time dependent. 
The phonon bath Hamiltonian has the usual form (fi, = 1) 



(3) 



where the dispersion relation cjq is specified below. The 
electron-phonon interaction has the linear coupling form^ji^ 



(4) 



q 1=1 



where Ni is the number of excess electrons in the i-th dot and 
= Aq e-'i-^' P,(q), with Ri = and R2 = d the dot 
position vectors, see Fig. |2] The dependence of the coupling 
constant Aq on the material parameters and on the wave vector 
q will be specified below. The dot form factor is 



P,(q) ^ j <f 



•ni{v) e 



(5) 



where ni (r) is the excess charge density in the i-th dot. With 
no significant loss of generality, we will assume that the form 
factor is identical for both dots and, therefore, drop the i in- 
dex hereafter. In the basis {|1), |2)}, after dropping irrelevant 
constant terms, the electron-phonon interaction simplifies to 



Hsb = A'*, 



(6) 



where 



and 



'S^-E^q + (7) 
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FIG. 2: Geometry of the double quantum dot charge qubit. 

with gq Aq P(q) (l — e"*'* '*). The phonons propagate in 
three dimensions, while the electrons are confined to the plane 
of the underlying two-dimensional electron gas (2DEG). No- 
tice that the electron-phonon coupling is not isotropic for the 
DQD (Fig. O: Phonons propagating along = and any 6 
do not cause any relaxation, while coupling is maximal along 
(f) — — TT /2 direction. We neglect any mismatch in phonon 
velocities at the GaAs/AlGaAs interface, where the 2DEG is 
located. 

We now proceed with the Born-Markov-Redfield 
treatment'^* '"* of this system. While the Born approximation 
is clearly justified for weak electron-phonon interaction, the 
Markov approximation requires, in addition, that the bath 
correlation time is the smallest time scale in the problem. 
These conditions are reasonably satisfied for lateral GaAs 
quantum dots, as we will argue below. 

Let us assume that the system and the phonon bath are dis- 
entangled ait — 0. Using Eqs. and (|5Jl, we can write 
the Redfield equation for the reduced density matrix p{t) of 
the DQDi^ 



p{t) = ~i [Hs{t),p{t)] + {[A(t)p(i), K] + H.c.} . (8) 

The first term on the right-hand side yields the Liouvillian 
evolution and the other terms yield the relaxation caused by 
the phonon bath. The auxiliary matrix A is defined as 



/•oo 

Jo 



(t) 



(9) 



where B{t) = Trfc{<I>(T)$(0)/(7?B)} is the bath corre- 
lation function, $(t) = e*^^'" $ e-*^«^, and /{Hb) = 
g-/3HB /Trb{e"''^«}, with P=l/T the inverse lattice tem- 
perature (kB = 1). 

Using Eq. Q in the definition of the bath correlation func- 
tion, we find that the latter can be expressed in the form 

poo 

B{t) = / dujiy{u;){e'^'^nB{Lo) + e-'^^[l + nB{uj)]}, 
Jo 

(10) 

where nsi'jj) is the Bose-Einstein distribution function and 



(11) 



is the spectral density of the phonon bath. 

We now specialize to linear, isotropic acoustic phonons: 
Cijq = s|q|, where s is the phonon velocity. Moreover, we 
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only consider coupling to longitudinal piezoelectric phonons, 
neglecting the deformation potential contribution. For bulk 
GaAs, this is justifiable at temperatures below approximately 
10 A'. Thus, 



2 2 

n\q\ ' 



(12) 



where (/ph is the piezoelectric constant in dimensionless form 
(5ph~0.05 for GaAs"'-'^) and fl is the unit cell volume. 

The excess charge distribution in the dots is assumed Gaus- 
sian: 



n{r) = S{z) 



1 



27ra- 



■ exp 



y 



2a2 



(13) 



This is certainly a good approximation for small dots with 
few electrons, but becomes less accurate for large dots. The 
resulting form factor reads 

P(q) =e-('?'+9?)'^'/2. (14) 

Note that this expression differs from that in Refs. "6 and|^ 
where a three-dimensional Gaussian charge density was as- 
sumed. 

Using Eqs. il2\ and il4\ . as well as the DQD geometry of 
Fig. 121 we get 



r/2 



d9 sin 9 exp 



Lud 



sin^e] (15) 



1 — Jn — sin 



or (ii) by changing the energy level difference e{t) keeping v 
constant (bias pulsing). Tunnel pulsing seems advantageous 
as it implies fewer decoherence channels and less leakage. 
However, a recent experiment used a bias pulsing scheme>ii 

Our system's Hilbert space is two-dimensional by construc- 
tion [see Eq. hence there is no leakage to states outside 
the computational basis. We can, therefore, use square pulses 
instead of smooth, adiabatic ones. This not only allows us to 
analytically solve for the time evolution of the reduced den- 
sity matrix, Eq. (|8}, but also renders our results applicable to 
both tunnel and bias pulsing. Indeed, in both regimes one has 
e{t) ~Q and v{t) — Vm for t > 0, taking that the pulse starts 
at t = 0. Let us assume that the excess electron is initially 
in the left dot: /Oii(O) — 1 and puiO) — 0. In this case, since 
the coefficients on the right-hand side of (jSji are all constants 
at t > 0, we can solve the Redfield equation exactly (see Ap- 
pendix|A|for details). As p{t) has only three real independent 
components, the solution is 



Pii(t) 
Repi2(t) 
Impi2(t) 



^ + ^\cos Lut + ^ sin ujt), (18) 
2 2 2lo 



-l(l-e-'^*)tanh^, 
e 2 smujt, 

ZUJ 



where 



(19) 
(20) 



It is instructive to inspect the asymptotic limits of this equa- 
tion. At low frequencies, 0) w lo^ /Qs^; thus, 
the phonon bath is superohmic. At high frequencies. 



{ \ ^ 9phS' ( d 



where 



(16) 
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(21) 
(22) 
(23) 
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dx X e ^ 
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Jo I -X 

a 



(17) 



Notice that the spectral function does not have the exponential 
decay familiar from the spin-boson model, but rather falls off 
much more slowly: v{uo — > oo) cx uo^^. This should be con- 
trasted with the phenomenological expressions used in Ref.0. 

The characteristic frequency of the maximum in v{ijj) is 
TjT^ = s/a. For typical experimental setups, a«50 nm while 
s w 5 X 10"^ m/s for GaAs, yielding Tc ~ 10 ps (r~^ « 65 ^eV). 
Thus, the Markovian approximation can be justified for time 
scales t > Tc and if all pulse operations are kept adiabatic on 
the scale of Tr- 



ill. DECAY OF CHARGE OSCILLATIONS 

One can operate this charge qubit in two different ways: (i) 
by pulsing the tunneling amplitude v{t) keeping e constant. 



Note that 71,2 <C v„i. We extract the customary energy and 
phase relaxation times, Ti and T2, by rotating to the energy 
eigenbasis {|— ), |+)}: 



(24) 



P-+it) 



^^ + Pii{t)+ilmpi2{t). (25) 



Then, the damping of the oscillations in the diagonal ma- 
trix elements is the signature of energy relaxation, while the 
phonon-induced decoherence is seen in the exponential decay 
of the off-diagonal elements. For the DQD, we find Ti—^^^ 
and T2 = 27^^ for the decoherence time. 

The quality factor of the charge oscillations in Eq. ( II8I 1 is 
Q = uj/ttji. Using Eqs. (12 1> . (I22> . and il5\ . we find that 



4 



4tanh(?;„/r) 



da; 



1- Jo 



(26) 



r 



where uja — s/2a. The Q-factor depends on the tunnehng 
amplitude u,„, lattice temperature T, dot radius a, and interdot 
distance d. 

Several experimental reaUzations of DQD systems recently 
appeared in the literatureiii*i2iiS^ In principle, all these se- 
tups could be driven by tunnel pulsing to manipulate charge 
and perform single-qubit operations. To understand how the 
Q-factor depends on the tunneling amplitude Vm in realis- 
tic conditions, let us consider the DQD setup of Jeong and 
coworkerspi^ In their device, each dot holds about 40 elec- 
trons and has a lithographic diameter of 180 nm. The ef- 
fective radius a is estimated to be around 60 nm based on 
the device electron density. Therefore, d/a « 3. The lattice 
base temperature is 15 mK. Introducing these parameters into 
Eq. i26\ . one can plot Q-factor as a function of u,,, or, equiv- 
alently, as a function of the period of the charge oscillations 
P = 2TT/ujfiiTT/vm- This is shown in Fig.|3l 

To stay in the tunnel regime Vm should be smaller than the 
mean level spacing of each QD, approximately 400 fieV in 
the experiment,'^ Therefore, in Fig.|3]we only show the curve 
for Vm up to 100 fj,eY. One has to recall that at these values 
the Markov approximation used in the Redfield formulation is 
not accurate (see end of Sec. HJ, and so our results are only 
an estimate for Q. For strong tunnehng amplitudes, when 
25 ^eV < Vm < 100 ^cV, the largest value we find for Q 
is close to 100. For weak tunneling with Vm < 25/ieV, the 
situation is more favorable and larger quality factors (thus rel- 
atively less decoherence) can be achieved. Nevertheless, the 
one-qubit operation time, which is proportional to the period. 
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FIG. 3: The charge oscillation Q-factor as a function of the tunnel- 
ing amplitude Vm (lower scale) and of the oscillation period P (upper 
scale) for a GaAs double quantum dot system. The lattice tempera- 
ture is 15 mK and the dot radius and interdot distance are 60 nm and 
180 nm, respectively. The inset shows the relation between P and Q 
at small tunneling amplitudes (large periods). 



grows linearly with Q in the region of Vm ^ 0, as shown in 
the inset of Fig. |3] Therefore, at a certain point other deco- 
herence mechanisms are going to impose an upper bound on 
Q. 

The minimum of Q in Fig.|3loccurs when Vm coincides with 
the frequency at which the phonon spectral density is maxi- 
mum. It corresponds to the energy splitting between bond- 
ing and anti-bonding states of the DQD, 2vm, being approx- 
imately equal to the frequency of the strongest phonon mode 

s/a: Vm ^ i^a- 

From Fig. [Sj it is evident that one can reach certain values 
for the Q-factor (say, Q = 100) at both weak {vm — 4.6 /leV 
~ 53 mK) and strong {vm — 93 /xeV ~ 1.1 K) tunnel- 
ing. However, these two regimes are not equally convenient. 
From Eq. \26\ . it is clear that the temperature dependence of 
the Q-factor is fully determined by the bonding-antibonding 
splitting energy 2t)„i: (5(T) = (5(0) tanh(t)„i/r). We no- 
tice that Q{T) « (3(0) if r < therefore, the (5-factor 
is less susceptible to temperature variations for strong tunnel- 
ing (Fig.l^. Another parameter that influences the (J-factor is 
the dot radius, which controls the frequency of the strongest 
phonon mode, s/a. In the strong tunneling regime (dashed 
curve in the inset to Fig. 13, one has to increase the QD size 
to improve the (J-factor This would reduce the energy level 
spacing, hence only moderate improvement in (J-factor is pos- 
sible. In contrast, in the weak tunneling regime (solid curve 
in the inset to Fig.|4} one has to reduce the QD size. This can 
lead to a significant (up to one order of magnitude) (J-factor 
improvement. 




0.4 0.6 
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FIG. 4: The charge oscillation (J-factor as a function of the lat- 
tice temperature. Inset: as a function of the dot radius for a fixed 
ratio d/a — 3. The solid (dashed) line corresponds to the weak 
Vm — 53 mK (strong v,n ~ 1.1 K) tunneling regime. Other parameter 
values are equal to those in Fig.|3| 
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IV. BIAS PULSING 



In a recent experimentpi Hayashi and coworkers stud- 
ied charge oscillations in a bias-pulsed DQD. In this regime 
the energy difference between the left and right-dot single- 
particle energy levels is a function of time: £{t) — £q u{t). A 
typical profile used for pulsing is 



If , i + WI2 
u{t) ^ 1 - ^ (^tanh 



t-Wl2\ 
tanh — j , (27) 



2r 



where W represents the pulse width and r controls the rise 
and drop times. During bias pulsing, the tunneling amplitude 
is kept constant. In Ref. 1 1 , the difference in energy levels 
was induced by applying a bias voltage between left and right 
leads (and not by gating the dots separately). For their setup, 
the maximum level splitting amplitude was Eq ~ 30 /xeV and 
T w 15 ps, corresponding to an effective ramping time of about 
100 ps.^" The tunneling amplitude was kept constant and es- 
timated as D « 5 /icV, which amounts to charge oscillations 
with period P w 1 ns. The lattice temperature was 20 mK. 
Each quantum dot contained about 25 electrons and the ef- 
fective dot radius is estimated to be around 50 nm based on 
the device electron density. From the electron micrograph of 
the device one finds d « 225 nm, hence d/a k, 4.5. When 
substituting these values into Eq. M6\ . one finds Q w 54. 

However, from the experimental data one observes Q w 
3. Low Q-factors were also obtained by Petta and cowork- 
ers in an experiment where coherent charge oscillations in a 
DQD were detected upon exciting the system with microwave 
radiationii^ Other mechanisms of decoherence do exist in 
these systems, such as background charge fluctuations^' and 
electromagnetic noise emerging from the gate voltages. Our 
results combined with the recent experiments indicate that 
these other mechanisms are more relevant than phonons. 

We now turn to yet another possible source of decoherence: 
Leakage to the leads when the pulse is on.^^ To illustrate this 
alternative source of damping of charge oscillations, we sim- 
ulate the bias-pulsing experiment of Ref. 1 1 by implementing 
a rate equation formalism similar to that used in Ref. 23. The 
formalism is based on a transport theory put forward for the 
strongly biased Mvmir^^ First, we find the stationary current 
/o through the DQD structure when the pulse is off (that is, 
the bias is applied) :^^ 



Lrr 



Li R 



(28) 



where e is the elementary charge. Ti^^r^ is the partial width 
of the energy level in the left (right) dot due to coupling to the 
left (right) lead (when the bias is applied); in the experiment, ' ' 
FL.i? ~ 30 /ieV. On the other hand, when the pulse is on, the 
stationary current is zero. We now apply the pulse e{t) and 
measure the current I{t). In the experiments, the level widths 
^L,R decrease upon biasing the system. To include that effect 
here, we also pulse them: Tj^it) =Jl + {^l ^ lL)u{t) and 
analogously for Tr, where ^l(r) is the residual leakage to the 
left (right) lead when the pulse is on. We use ^l.r = 0.3 /xeV, 
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FIG. 5: (a) The response current /rcsp(i)/e in ns^^ as a function 
of time for a pulse with = 4 ns and t = 30 ps. (b) Number 
of electrons transferee! between left and right leads, as defined in 
Eq. |29j, as a function of the pulse width W . 

even though the real leakage in the experiment was likely 
much smaller To obtain the response current one subtracts 
the stationary component: Ircs-p{t)—I{t) — Iou{t). 

FigurelSja) shows the response current for a pulse of width 
= 4 ns and r = 30 ps. The latter is approximately twice as 
large as in the experiment and is chosen to enhance the effect. 
In Ref. 11, pulses were applied at a frequency / — 100 MHz. 
The average number of electrons transfered from the left to the 
right lead per cycle minus that in the stationary regime isS 



1// 



dt Ircsp{t)/e. 



(29) 



(In the simulations there is no need to apply a sequence of 
pulses.) Notice that n oscillates as a function of the pulse 
width W [see Fig.|5jb)] as observed in the experiment. Two 
main conclusions can be drawn from our simulation. First, the 
larger r, the smaller the visibility of the charge oscillations.^— 
Second, the larger the leakage rates jl.r^ when the pulse is 
on, the stronger the damping of the oscillations. While the 
damping due to leakage is presumably too weak an effect to 
discern in the data presented in Ref. 1 1 , the loss of visibility 
due to finite r is likely one of the causes of the small amplitude 
seen experimentally. 



V. CONCLUSIONS 

The main conclusion of the paper is that, under realistic 
conditions, phonon decoherence is one to two orders of mag- 
nitude weaker than expected. ''■^•^ The analytical expression for 
the Q-factor given in Eq. ( I26> was found using an expression 
for the phonon spectral density, Eq. ( I15t . which takes into ac- 
count important information concerning the geometry of the 
double quantum dot system. In a previous work^ an approxi- 
mate, phenomenological expression, i^{uj) oc w exp(— w/wc), 
was utilized in the treatment of charge qubits. There is a 
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striking difference between these two expressions in both the 
high- and low-frequency Hmits. Moreover, an arbitrary cou- 
pHng constant was adopted in Ref. 7 to model the electron- 
phonon interaction while our treatment uses a value known 
to describe the most relevant phonon coupling in GaAs. On 
the other hand, other previous work^ assumed a spherically 
symmetric excess charge distribution in the dot while we have 
assumed a two-dimensional pancake form. These differences 
account for most of the discrepancy between the present and 
previous results. 

Based on these findings we conclude that phonon decoher- 
ence is too weak to explain the damping of the charge oscil- 
lations seen in recent experiments." '^ Charge leakage to the 
leads during bias pulsing is an additional source of damping, 
as shown in Fig.|5jb); however, for realistic parameters ; it 
turns out to be a weak effect as well. Hence, other decoher- 
ence mechanisms, such as background charge fluctuations or 
noise in the gate voltages, play the dominant roleiSi 

There are two distinct ways to operate a double quantum 
dot charge qubit: (i) by tunnel pulsing or (ii) by bias pulsing. 
Tunnel pulsing seems advantageous due to the smaller num- 
ber of possible decoherence channels. In addition, the bias 
pulsing scheme, in contrast to tunnel pulsing, introduces sig- 
nificant loss of visibility in the charge oscillations. 

In this work we did not attempt to study leakage or loss 
of fidelity due to non-adiabatic pulsing, which are both im- 
portant issues for spin-hased quantum dot qubits.^^ Moreover, 
we have not attempted to go beyond the Markov approxima- 
tion when deriving an equation of motion for the reduced den- 
sity matrix. Both of these restrictions in our treatment impose 
some limitations on the accuracy of our results, especially for 
large tunneling amplitudes. 

Finally, it is worth mentioning that some extra insight 
would be gained by measuring the Q-factor as a function of 
the tunneling amplitude u„, experimentally. Such a measure- 
ment would allow one to map the spectral density of the bo- 
son modes responsible for the decoherence. This would pro- 
vide very valuable information about the leading decoherence 
mechanisms in double quantum dot systems. 
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FAFERJ, and PRONEX. 



APPENDIX A: DERIVATION OF EQS. gS), {T5), AND ^ 

For t > 0, Eq. (|3 is time-independent: Hs — WmC^- Since 
the K matrix is also time-independent [Eq. (0], the matrix A 
defined by Eq. (|9} is time-independent as well. After some 
straightforward operator algebra, we find that 

A = i / dT B{t) e-"^rr.o^ ^irv^<y. (Ai) 

2 Jo 

1 [°° 

= - / dr B{t) [oz cos.{2vmT) - Gy sin(2wmr)] . (A2) 
^ Jo 

One can rewrite Eq. ( IA2> as follows 

^ = ^(71 + »73)ctz - ^(72 + «74)ctj;, (A3) 
where {7i}'s are real coefficients: 

/ 71 + ^73 1 _ (°° , 0/ N / cos(2v„r) 1 

The density matrix p{t) is a 2 x 2 Hermitian matrix with 
unit trace. Hence, it has three real independent components 
and can be written as follows: 

1 , 1, 

/? = - + fjj; Re P12 - crylmpi2 + (pn - -). (A5) 

Let us substitute Eqs. ( IA5t and iA3\ into the Redfield equation 
[Eq. (|8}] and use that Hs = VmCTx and K = \crz- A simple 
algebraic manipulation leads to three differential equations, 

74 

Re/9i2 = -7iRepi2 + Y, (A6) 

Pll = -2wmlmpi2, (A7) 

Im/?i2 = (2t;„i + 72)(pii - -) - 7ilmpi2. (A8) 

The initial conditions are pii{0) — 1 and pi2(0) ~ 0. Eq. ( IA6> 
decouples from Eqs. ( IA7> and ( IA8t . Its solution is given by 
Eq. (I19> . where we used the following identity: 74/71 — 
— tanh(w,„/T). Eqs. ( IA7t and ( IA8> form a closed system. 
Their solution is given by Eqs. ( I18> and ( I20> . 

The coefficients 71 and 72 [Eqs. M2\ and (I23> . respectively] 
are calculated using Eqs. (IA4> and jlO> . 
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